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The plastic flow of a polycrystal is analyzed assuming grains as fine that the rate limiting pro- 
cess is grain boundary sliding, and grains readily accommodate their shapes by slip to preserve 
spatial continuity. It is shown that thinking of a polycrystal with randomly oriented grains as an 
homogeneous and isotropic continuum when dealing with it as a dynamical medium, even in a scale 
much larger than the grain size, leads to gross errors. The polyhedral nature of grains influences 
the plastic flow in a radical manner, as the relative velocity of adjacent grains is constrained to the 
common boundary plane, and only the in-plane shear stress contributes to their relative motion. 
This constriction determines that V ■ v 7^ 0, where v is the velocity field of the material medium, and 
plastic deformation necessarily involves grain volume variations, which can only be elastic. As this 
has a limit, fracture follows as a necessary step of plastic flow. A theoretical approach to plastic flow 
is developed, and emphasis is done in superplastic deformation, from zero to fracture strain. The 
theory allows to quantitatively explain the observed features of superplastic materials and responds 
to many open questions in the field. 

PACS numbers: 62.20.fg, 46.35. +z, 62.20.mj, 62.20.mm 



I. INTRODUCTION 

Superplasticity is an anomalously high tensile ductil- 
ity occurring in a growing class of fine grained polycrys- 
talline solids. At rather precise conditions of composi- 
tion, grain refinement, temperature and strain rate, they 
can undergo uniform neckless plastic elongations that 
may reach thousands of percent prior to fracture. These 
large strains occur with no significant change in the grain 
shapes or sizes, nor in any other fine attribute [TJ[2]. As 
well, different bodies in the superplastic regime of de- 
formation become prone to diffuse one into the other 
when their surfaces make contact (diffusion bonding). 
Among the many singular features of the phenomenon 
is the fairly narrow range of temperatures over which it 
happens. For instance, tensile tests for the superplas- 
tic titanium alloy Ti-6A1-4V give maximum ductility at 
temperatures close to 1153 K, and elongation to frac- 
ture rapidly decays to half the optimal figure within a 
neighborhood not wider than 200 K around [3]. These 
temperatures are well below the melting point (1930 K), 
so we are talking about a genuine solid state effect. 

Beside the scientific problem it poses, superplasticity 
has acquired considerable technical importance. After 
the pioneer work of Backofen on the superplastic form- 
ing of Zn-Al sheets by a simple air pressure operation 0], 
as in the fabrication of plastic bottles, and the discovery 
of a number of high strength superplastic alloys, the sub- 
ject moved into the scope of engineering materials for ad- 



vanced applications. Since then, the list of high strength 
superplastic alloys, intermetallic compounds and ceram- 
ics has increased notably [2], together with the technical 
knowledge for producing superplastic structures. A num- 
ber of present day aluminium superplastic alloys, like Al- 
7475, Al-8090 and Supral 220, or titanium Ti-6A1-4V and 
SP700, have outstanding importance as aeronautic ma- 
terials [5 . Superplastic forming combined with diffusion 
bonding is now a standard technique for producing both 
engine and airplane structural parts. The development of 
light ductile materials of high strength at elevated tem- 
peratures has a principal role in the race for more efficient 
aeroengines. Superplastic titanium and nickel aluminides 
are drawing great attention in this respect [21 [5] . Having 
half the density of the current materials, TiAl-based al- 
loys are promising candidates to replace the nickel-based 
superalloys in combustion chambers and turbine blades. 

After a large plastic deformation with no significant 
microstructural changes, and resembling a steady flow, 
superplastic solids collapse at a rather well-defined strain, 
which depends on temperature and strain rate. The reg- 
ular behaviour of superplastic elongation to fracture sug- 
gests that the phenomenon is governed by a specific vari- 
able, not identified so far, which keeps account for the 
total strain and produces failure when reaching a pre- 
cise value. Certainly, the development of cavities may 
contribute to precipitate collapse [5J [7] but seems not 
to be the main cause. Maximum elongation to fracture 
strongly depends on the strain rate, while the evolution of 
cavities, when occurring, has been observed to be almost 
independent of that variable [HIS]- Moreover, titanium 
superplastic alloys do not cavitate significantly [3l [10] . 
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For decades, people have been interpreting their ob- 
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servations on the basis of constitutive equations of the 
form a n — e m /A, relating the uniaxial applied stress 
tr with the superplastic strain rate e. The coefficient 
A depends on the microstructure and follows an Ar- 
rhenius dependence with temperature. A variety of 
theories give this kind of equation with constant m, 
[IH H3 HH [T3 HEl H71 HE] - Nevertheless, experiments 
show that log a and loge are never linearly related, as 
implied by the power law, but the curves display a sig- 
moidal shape, which is a characteristic feature of the su- 
perplastic deformation, and the value of the strain rate 
sensitivity m that fits the data near the inflexion point is 
usually a non-integral number which vary strongly with 
strain rate and strain. In many superplastic materials 
the large ductility is observed even for log e far from the 
inflexion point, where the slope m has changed, and the 
strain rate at the inflexion point shifts considerably with 
strain [19] . 



Despite the success of some previous work in accurately 
predicting the temperature dependent stress vs. strain 
rate curves at zero or small strain for a number of su- 
perplastic solids [2H1 HU [22 , it rests to explain why are 
they superplastic. Knowing the mechanisms that make 
the grains to slide past each other and the successful 
prediction of plastic properties at small strains are not 
enough to answer this question, and it is necessary to in- 
vestigate how grains flow over long paths when the solid 
undergoes a finite deformation. The literature shows sev- 
eral different theoretical attempts to explain superplastic 
phenomena, but no one addresses the point of the modi- 
fications occurring inside the material when strained up 
to fracture ETJ E2J E2J El ESI El EH ESI 123 HI1 HI] , and 
the dependence of the stress vs. strain rate curves with 
strain. 



We put forward here a general theoretical approach 
for the plastic deformation of fine grained solids up to 
fracture strain. We assume the grain boundary sliding 
mechanisms of Refs. |2"T1 and 12^1 but refine the study of the 
grain collective motions in the spirit of Ref. 23 . The latter 
is also modified in order to interpret better the physics 
of the sliding grains. As a general conclusion, our results 
prove that the deformation and fracture properties of su- 
perplastic solids are governed by general laws amenable 
of analytical mathematical description, simple and accu- 
rate enough to constitute a practical tool for interpreting 
material tests. Though the assorted conglomerate of pos- 
sible structural defects may affect in specific situations, 
the main mechanical behaviour of solids is determined by 
a few parameters associated to the most basic structural 
properties. 



II. STRESS TENSOR AND PLASTIC STRAIN 
RATES 

A. General considerations 

At a scale much larger than the grain size, polycrys- 
tallinc matter lacks symmetry constrictions and period- 
icity, and displays same average packing and properties 
in all directions, and over its whole extention. However, 
despite this, assimilating a polycrystal to an homoge- 
neous and isotropic continuum may lead to gross errors, 
no matter the scale, when dealing with it as a dynamical 
medium. It has been proven before |23 , and is revised 
in this section, that the faceted nature of the structural 
constituents of a polycrystal, the grains, influences the 
macroscopic plastic flow in a radical manner. 

To explain briefly how this does happen, assume that 
plastic flow takes place by just grain sliding, and that 
grains readily accomodate their shapes by slip, so the 
rate limiting process is the former. Sliding of two adja- 
cent grains occurs when the shear stress in the plane of 
the common boundary exceeds a critical value t c . How- 
ever, no matter how big the externally applied forces may 
be, shear stresses vanish in planes whose normals are in 
the principal directions of the stress tensor. By conti- 
nuity, shear stresses in planes whose normals are inside 
well-defined finite solid angles around the principal di- 
rections are smaller than r c . Hence, many grains are im- 
peded to slide simply because of the orientation of their 
surfaces [52j. By this effect, the material medium ceases 
to be dynamically isotropic, and cones are created by 
the very external forces within which no displacement is 
allowed. These forbidden cones destroy the kinematic 
balance and V • v ^ 0, where v represents the velocity 
field of the flowing medium. Hence, mass density is not 
locally conserved in the overall plastic flow. The conse- 
quent internal pressure build up helps deformation and, 
at a critical strain, the material becomes mechanically 
unstable, undergoing brittle or ductile fracture, depend- 
ing on the precise physical circumstances. The internal 
pressure p, which operates in the local microscopic scale, 
keeps account of cumulated strain and triggers fracture 
at a precise value of it. 



B. Relative velocities between adjacent grains and 
strain rates 

The proper configuration variables in the macroscopic 
scale are the components £y = x, y, z) of the strain 
tensor, which, together with their time derivatives iij, 
give a complete picture of the dynamical state of the 
system. However, the evolution of in time is governed 
by forces excerted between each pair of adjacent grains, 
and hence the dynamical analysis demands to go first to 
the grain level. As the objective is to return to the larger 
scale, the first task is to establish the connection between 
the two descriptions. 
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Given the statistics of spatial orientations of the grain 
boundaries, which may be assumed isotropic from the 
beginning, one can derive the relation between and the 
set of the relative velocities Aw between pairs of adjacent 
grains. To show this, temporarily place the origin O of 
the coordinate system at a point of a grain surface. Then 
call v the velocity of a point P in the positive i-coordinate 
axis, which is in another grain boundary distant Xi = nd 
from O, where n is a large natural number and d the mean 
grain size. (In what follows it will be denoted either i,j = 
x, y, z or i, j — 1,2, 3, with x x = x, x-i = y, X3 = z). 

As there are n contiguous grains between points O and 
P, v is the sum 



77=^Aw(p) + ^u l (p)e l , (1) 
p— 1 p— 1 

where p correlatively numbers the adjacent grains along 
the i-axis and Av(p) is the velocity of grain p relative to 
grain p — 1 at the interface between them. The second 
term in the right hand side of Eq. accounts for the 
assumed deformability of the grains, ej is the unit vector 
along the i-axis and Ui (p) stands for the expansion rate of 
grain p along the i-direction. However, we are consider- 
ing plastic deformation processes in which grains remain 
equiaxed in the mean and do not vary significantly their 
average size. Therefore, on this basis, as n is large we 
can set 



implicit in the demonstration given in Ref. but that 
they conserve mean size and shape in the plastic flow). 

If the material is homogeneous and isotropic in the 
macroscopic scale, the n plane grain interfaces intersect- 
ing OP are oriented in any direction with the same prob- 
ability. Therefore, one may interpret the sum in Eq. {§) 
as an average which, as n is large, runs over all possi- 
ble orientations of the grain boundaries intersecting OP. 
However, vectors Av represent relative velocities between 
consecutive grains ordered along the positive z-coordinate 
axis. Thus, to take this into account, the average must be 
performed over all orientations in the hemispace Xi > 0. 
Equivalently, denoting h the unit vector normal to the 
interface between two grains, the sum in Eq. g can be 
replaced by an integral over all directions h satisfying the 
condition n ■ e, > 0. Therefore, we can write 

iij = j{Avj)i, (6) 

where symbol (...)» means the average over all orienta- 
tions in the semispace x% > 0. If the material is not 
isotropic in the larger scale, for instance because the 
grains are not equiaxed, the average ^ should incor- 
porate the proper weigth factors. 

C. Relative velocities between adjacent grains and 
stresses 



n 

Y,Mp)^ = 0, (2) 
P =i 

which means that, statistically speaking, the grains nei- 
ther stretch nor shrink. Hence 



*=x;mp). (3) 

If the solid is being deformed, the components of the 
displacement of P in a time St are 

ViSt = euXiSt, v j St=e ji x i St (j^i). (4) 

Substituting Xi — nd and Eq. ^ in Eqs. Q, it is ob- 
tained that 



p=l 

where the i dependence in the right hand side is implicit 
in the direction along which the relative velocities are 
sampled. (Notice that the present arguments leading to 
Eq. ([5]) do not assume that the grains are rigid, as seems 



The stress tensor 

/ a x \ 
K) =0^0, (7) 

V <7 Z ) 

is diagonal in the main frame of reference (xyz), with 
a - , (i = x,y,z) being the principal stresses. To properly 
express the physics at the grain scale, define a local frame 
of reference {x'y 1 ' z') whose x'y' plane is in the common 
grain boundary plane of two adjacent grains (Fig.[T]). In 
the local reference system the stress tensor reads 

{a ef ) = R(a i:j )R T , (8) 

where R is the rotation matrix that makes the transfor- 
mation from the (xyz) to the (x'y' z') frames, and R T is 
the transposed of R. 

As pointed out in Ref. 1211 superplastic grain boundary 
sliding seems closely related with elastic instabilities oc- 
curring at interfaces in solids when subjected to overcrit- 
ical stresses gll (111 H3 HH1 HH EH1 EH E21 E3] • Monte 
Carlo simulation and analytical treatment show that, in 
general, planar structures inside the bulk of crystalline 
matter become rough when subjected to high enough in- 
plane shear forces [23]. After the occurrence of the elastic 
instability the locally induced normal stress field drives 
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grain boundary 



o 



Avi 



= (*3. 
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■(x'y') 



FIG. 1: Local reference system (x'y'z'), with the z' axis 
normal to the plane of the common boundary of two adjacent 
grains. The axes of the (xyz) frame of reference are in the 
principal directions of the stress tensor. 



an active atomic transport between the interface and the 
two adjacent crystals. This bulk effect is closely related 
with the Asaro-Tiller-Grinfeld instability, i. e. stress in- 
duced roughening of solid surfaces [25] [26j [27] , and with 
the spontaneous roughening of thin films produced by 
strong in-plane stresses arising from lattice mismatch 
with the substrate [Ml HO US HQ GUI EB] . 

It has been shown that the active atomic transport be- 
tween the buckled grain boundary and the two adjacent 
grains follows closed paths, which makes the two grain 
surfaces to slide [21]. The relative velocity between the 
two sliding grains turns out to be proportional to the 
in-plane components, a x ' Z > and cr y ' z >, of the shear stress 
operating in the shared grain boundary plane, provided 
that the total in-plane shear stress 



>y'z> 



(9) 



be higher than the threshold stress r c for the elastic in- 
stability of the grain boundary [21, 22, 23 . Hence 



Aw 



0, 



T (x'y>) 



> T r . I = 



x',y', 



otherwise, 



Av z > = 0, 

(10) 

where Q is a coefficient depending on the normal stresses 
Oi'i' which will be discussed in detail in the next section. 
The factor (1 — t c /t z >) ensures that Av = for t z > = t c . 

The coefficient Q must not depend on the orientation 
of the local frame of reference (x'y' z'), therefore its de- 
pendence on the normal stresses is only via the invariant 



. (&x'x' 



;(a x + a y + a z ). (11) 



FIG. 2: The solid line represents the linear relation assumed 
here between |Av| and the total in-plane shear stress. The 
dotted line represents the model of Ref . [23] . 



This is also a necessary and sufficient condition for con- 
sidering Q(p) as independent of the shear stresses for any 
orientation of the grain boundary. But for the critical 
condition incorporated in Eqs. ( 10 ), Q may be seen as an 



inverse grain boundary viscosity coefficient [T3 [Ml US . 
Notice that Avy in Eqs. (10 1 is related with an irre- 



versible shift between the two grains, and hence cannot 
describe the variation of an elastic degree of freedom. In 
what follows, elastic and plastic strains will be indepen- 
dent mechanical variables (which nevertheless do inter- 
act). Grain clastic distortions arc always implicit in the 
subsequent developments, but plastic strains are the dy- 
namical variables because grains are the flowing entities. 
Elastic strains of the grains will remain implicit in the 
stresses, which are linearly related with elastic strains 
through the equations of the theory of elasticity. 

To put Eqs. ( 10 1 in a more practical way, we must 
perform the matrix transformation ^ , choose the proper 
components of (cr^ji) to insert them in Eqs. (10 1 and 
write Av in the local (x'y'z') reference system, and then 
make the inverse vector transformation to express Av in 
the main (xyz) coordinate system. Hence, for t z i > t c 
the components Avi> are 



Av v = Q[R(a ij )R T ]., z , ( 1 
Av z > = 0, 



r (x'y') 



i' = x',y', 



(12) 

where [■■•] i / z / selects the i'z' component of the tensor. 
Wh en expressed in the (xyz) frame, the relative velocities 
(12 1 are (Avi) = R T (Avi') or, more explicitly, 



Avi = ^ Ri'iAv^, 



x,y,z. 



(13) 



The local (x'y'z') frame can be arbitrarily rotated 
around the z' axis without breaking the condition that 
z' be normal to the grain boundary plane. Then, the x' 
axis can always be chosen on the xy plane, as in Fig. [3] 
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Calling 4> the Euler angle between the x' and x axes, and 
8 the other Euler angle, going from the z to the z' axis, 
the rotation matrix reads 



R{6, 



cos <p 
- sin cj> cos ( 
sin d> sin 8 



sin cf) 
cos (f> cos sin 8 
— cos </> sin 8 cos 



(14) 



-((9, <j>) = [{a x - cr y ) 2 sin 2 0cos 2 </> 



sin 

1/2 ( 18 ) 

+ (a x sin 2 </> + <t v cos 2 - u z ) 2 sin 2 cos 2 0] , 
otherwise Au = 0. 



D. Strain rates and stresses 




To simplify the equations, let us particularize to the 
important case of uniaxial external stress a z . Far from 
the sample surfaces, the material has cylindrical symme- 
try and a x = a y , which will be assumed not to vanish for 
the sake of generality. Eqs. (16 1 and (18) greatly reduce 
and give 



Av = Q (p z — <r x ) sin 8 cos ( 
xll 



V {.Ox = <Jy), 



(19) 



(a z — <j x ) sin 8 cos ( 
where the unit vector v in the {x, y, z) system is 



FIG. 3: Angles between the axes of the local reference system 
(x'y'z') and the principal directions. The x' axis is in the xy 
plane. 



Transforming (cr,j) with the matrix (14 1 one obtains, 
in particular, that 



&x'y' — Wx ~ &y) sin (/jcos </> sin ( 



>y'z- 



-(a x sin 2 cj) + (Ty cos 2 <fi — a z ) sin 8 cos 8, 



(15) 



and, making the inverse vector transformation (13 1, the 



components of the relative velocity between grains takes 
the explicit form 



v = (— sin (p cos 8, cos cj) cos 8, sin ( 



(20) 



The choice of the x' axis in the xy plane proves to be 
particularly convenient because, by the cylindrical sym- 
metry, any direction in the xy plane is a principal direc- 
tion, and then a x i z i = for any 8 and 4>. The condition 



( 17 1 on the minimal shear stress r c becomes 



\(<t z - a x )sm8cos8\ > t c , (<j x = <j y ) 7 (21) 
which determines a critical angle 8 C1 given by 



sin(26» c ) 



2t> 



(22) 



Au x = Q sin 8[(a x — a y ) sin cos 2 4> + (a x sin 2 
+ (T y cos 2 cj) ~ Oz) sin </>cos 2 8] ( 1 TK~T 



Av y = Qsm8[(a x — a y ) sin 2 </>cos 4> — (a x sin 2 

+ (T y COS 2 (j) — <T Z ) COS (j) cos 2 8] ( 1 JW~, 



Av z = —Q((T X sin 2 



(TyCOS <p — a z ) sin 8cos8 1 



r(0,0) 



Eqs. ( 16 I hold when 



t(0,<£) > T c , 



with 



(16) 



(17) 



such that Av = for | sin26'| < sin2# c . Therefore, under 
the restriction that a x — <r y , plastic flow is not allowed for 
< 8 < 8 C , vr/2 - 8 C < 8 < tt/2 + 8 C and tt - 8 C < 8 < n, 
as shown in Fig. |4ja). 

To calculate the averages (Avj)i, appearing in Eq. pi, 
we must integrate over all directions (8, cf>) such that the 
unit vector h{8 1 </>), normal to the grain boundary plane, 
has positive projection on the Vector h is the 

unit vector k' along the z' axis of the local coordinate 
system, that is, 



h = (sine/) sin 8, — cos <^>sin 8, cos 8). (23) 



In general, 



{Av) t = —j dcj)d8 sin 8 Av(8,(t)), (24) 
27r Jft 




(a) 



(b) 



Eqs. (261 and ( 27 I hold for o z -o x > 0. When a z -a x < 



they must be modified reversing the sign of the terms 
multiplied by r c . 

By Eq. (22 1, a z — a x can be written in terms of just 
9 C and the constant t c . Replacing in Eqs. (26) and d27|, 



the strain rates turn out to be entirely determined by 
9 C , which becomes a natural and convenient auxiliary 
variable which varies in the range 





< C < it/A. 



(28) 



(c) 



(d) 



For finite r c , 9 C — corresponds to the asymptotic situ- 
ation in which the external forces applied on the sample 
are very strong. For 9 C — 7r/4, the forbidden directions 
extend over the whole sphere and no grain is allowed to 
slide. 

As the relative velocities, the strains are purely 
plastic, but elastic strains are not disregarded and re- 
main implicit in the stresses, which are linearly related 
with elastic strains through the equations of the theory 
of elasticity. 



FIG. 4: (a) Grain boundaries whose normals are in the spher- 
ical shells are the only ones allowed to slide, if the stress 
field has cylindrical symmetry with respect to the z axis. 
Grain boundaries whose normals make angles < 8 < 8 C , 
n/2 — 8 C < 9 < n/2 + 8 C and n — 9 C < 8 < n with the z axis 
do not contribute to plas tic flow, (b), (c) and (d) Integration 
domains for integral ( 24 1. 



III. 



DILATION RATE AND INTERNAL 
PRESSURE BUILD UP 



From Eqs. ( 26 I and ( 27 1 , the dilation rate V/V = e x 



turns out to be 



and the integration domains Di can be determined by 
examining the geometry of the {xyz) — * {x'y 1 ' z') trans- 
formation. They are (Fig. |4jb),(c) and (d)) 



D x :9e [9 C , tt/2 - C ] U [tt/2 + 9 c ,tt- 9 c ], 
<t>e [0,7r], 

D y : 9 e [9 C , tt/2 - 9 C ] U [tt/2 + 6 C , it - C ], 

[7r/2,37r/2], 
D z :9e [d c ,n/2-e c ], 
<be [0,24 



(25) 



Integrating and combining with Eq. ([6| one obtains 



Q 

id 



[(a z - o x ) cos(20 c ) - t c (tt - 40 c )] , (26) 



and 



V 
V 



sin(4^ c ) 



1 - cos(26> c ) - 



A9 r . 



- 29 c sin(26» c ) 



(29) 



where V stands for the volume of the sample. 



Eq. ( 29 1 discloses the amazing feature that the dilation 



rate does not vanish in general if the threshold stress for 
grain boundary sliding r c =/= 0. If r c = Eqs. (26 1 and 
(27 1 reduce to 



Q 



<y z — o x 



(30) 



-^(o- 2 -cr x ) (r c = 0), 



V/V — and the material flows as a viscous liquid. 

The dilation rate (29) vanishes only in the extreme 
cases 9 C = and 9 C = tt/4, and for a z — a x > is neg- 
ative for any other value in between. This means that 
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plastic stretching involves a spontaneous compression of 
the grains, which can only be elastic, and demands the 
ability of the material to support some clastic volume 
reduction. As this ability is essentially limited, plastic 
strain is limited as well, and a finite strain to failure is 



implicit in Eq. ( 29 1 . The single-crystal elastic compress- 



ibility will determine how much the polycrystal could be 
plastically elongated. 

Another general consequence of Eq. p9l is that, in 



general, the transversal stresses a x = a y do not vanish in 
the plastic deformation, even when the applied external 
force is strictly uniaxial, and the material is isotropic in 
a scale much larger than the grain size. This may seem 
rather odd if one has in mind a static, or quasistatic, pic- 
ture of the deformation process. Such a picture is, how- 
ever, misleading. Plastic flow is an essentially dynamic 
process. 



It must be noted that Eq. (29) deals with volume 



changes, which are necessarily elastic, but has nothing 
to do with grain shape variations. Our primary hypoth- 
esis is that grains are plastic, and they readily accom- 
modate their shapes to preserve matter continuity, but 
grain boundary sliding is the rate limiting process. This 
is consistent with accomodation by slip, since crystal de- 
formation by slip does not demand much additional stress 
once the critical resolved shear stress (CRSS) has been 
exceeded. Then, if t c is greater than the CRSS of the ma- 
terial then our hypothesis about grain plasticity should 
hold. We are not very specific about the accommoda- 
tion mechanisms because grain plastic deformation is not 
the rate limiting process, diffusional creep or dislocation 
creep may contribute, but we assume that the main ac- 
comodation mechanism is by slip processes. 

Settling <j x — cty — for e zz = is a natural initial 
condition. However, it will be shown that, as the plas- 
tic stretching along the z axis increases with time, the 
magnitude of the transversal stresses a x — a y monoton- 
ically increases as well, with negative sign. Hence, the 
plastic elongation of the sample produces a spontaneous 
compression in the transversal plane. Neck formation 
prior to ductile fracture of polycrystalline bars can be 
attributed to this circumstance. 



The dilation rate ( 29 ) can be understood in a qualita- 



tive manner by a simple geometric argument. The rela- 
tive speed | Aw | between adjacent grains depends only on 
9 and is symmetric with respect to 9 — 7r/4 in the top 
hemisphere, and with respect to 9 = Stt/A in the bottom 
one (Eq. (19l). Fig. [5] depicts the solid angles A, B and 
C in which Av = 0. Grain boundaries whose normals h 
are in B correspond to banned relative velocities Av in A 
or C, and vice-versa. As the solid angle B is larger than 
the sum of A and C, there are more supressed vectors 
Av in A and C than in B. Therefore, grain displace- 
ments in the z direction, or close to it, are less favoured 
than those occuring in directions near the xy plane. If 
no dilation takes place when 9 C — (i. e. r c = 0), some 
volume change should occur for r c ^ 0. The same kind 
of argument serves to show that no dilation does occur in 




FIG. 5: The forbidden directions for Av. Vectors Av in the 
solid angles A and C correspond to grain boundaries with n 
in the solid angle B, and vice-versa. Then, there are more 
vectors supressed in A and C than in B. 



two dimensions, and the area is conserved in the plastic 
distortion of a strictly two-dimensional layer. 

We will examine the case of the uniaxial plastic strain 
of a bulk sample along the z axis at a strain rate e = e zz . 
The only externally applied stress is a = a z . We will 
denote also p — —{cr x + a y + cr z )/3 for the pressure. The 
transversal stresses a x = a y are assumed finite to give a 
chance for internal stresses to develop. Of course, cylin- 
drical symmetry demands that the point into considera- 
tion be far enough from the surfaces. With this notation 

1 3 

o x = o y = + 3p), cr 2 - a x = -{a + p). (31) 



Recalling Hooke's law 

AV 

~V~ 



B 



P, 



(32) 



where B is the bulk elastic modulus, the dilation rate 
reads 



V 
V 



1 dp 
Bdi' 



(33) 



Replacing the identity edp/de = dp/dt and Eq. (29 1 in 
Eq. (f33j), 



dp _ Bt c 
de 2de 



29 A 1 



1 - cos(26» c ) 
sin(26» c ) 
2 

7rsin(2# c ) 



2cos(26» c ) 



(34) 



where use was made of Eq. ( 22 1 to eliminate a z — a x and 



put the right hand side in terms of just 9 C . The factor s 
is 



8 



if o z 
if a z 



<J X > 
cr x < 0. 



(35) 



For traction (s = 1), the right hand side of Eq. (34) 



is positive for < 8 C < 7r/4, and hence the pressure p 
increases monotonically with strain. 

Elastic degrees of freedom, implicit in V/V, are 
brought to the fore once more by Eq. (29). Eq. (33) 



expresses the volume variation (the trace of the elastic 
strain rate tensor) in terms of the pressure p (essentially 
the trace of the stress tensor) . Hence it allows to return 
to the scheme in which the elastic distortions are kept 
implicit in the stresses. The variation of p in the plastic 
deformation indicates that elastic and plastic degrees of 
freedom are coupled. 

Written in terms of just the critical angle 9 C , Eq. (26) 
for the strain rate becomes 



sQ 



2d 



cot(20 c ) + 26 c 



(36) 



stressed polycrystaline material must distinguish inter- 
granular and crystalline matter because they have dif- 
ferent mechanical properties [5T] [52]. A shear stress 
greater than a critical value applied in the plane of a grain 
boundary should buckle it, producing either a bound- 
ary corrugation or the formation of a periodic series of 
trenches. This localized deformation induces a periodic 
normal stress field that alternates compression and trac- 
tion on the adjacent crystal surfaces, as shown schemati- 
cally in Fig.[6ja). Grain boundaries are efficient sinks and 
sources for vacancies and the periodic stress field induced 
by the buckled boundary yields a periodic variation of the 
equilibrium value for the concentration of crystal vacan- 
cies. The grain boundary then evaporates and condenses 
the point defects in alternate sectors, producing streams 
of the defects in closed loops that cross the boundary 
and involve the two adjacent crystals (Fig. [6^b)). These 
closed loops perform as the driving pulleys of a conveyor 
belt and the phenomenon provides either an accommoda- 
tion mechanism and the driving force for crystal sliding 



with 



- arcsm 
2 



4r r 



3 k 



■p\ 



(37) 



In what follows we will consider just traction, and s will 
be omitted. 

At the start of the deformation process e = and 
&x = &y = 0. Then the initial condition written in the 
new notation is 



for 



0. 



(38) 



IV. THE COEFFICIENT Q 



The coefficient Q, governing grain boundary sliding, 
was studied in detail in Refs. [2T] and [55] . It was shown 
there that the stress dependence of Q is only through the 
pressure p, and varies with temperature and grain size. 
It obeys the rather simple law 



exp 



e + Q*P 



(39) 



where kg is the Boltzmann constant, T the temperature, 
the coefficient Co depends only on the grain size d, the 
constant eo is the energy necessary for evaporating a crys- 
tal vacancy from the grain boundary, and f2* measures 
the sensitivity of this energy to stress. 

The derivation of expression (39 1 has been already 



published [5TJ [55] , but a brief explanation of the physics 
involved is in order here. The mechanical analysis of a 



Traction Compression 



Grain boundary 




Compression Traction 
(a) 

Diffusion 



Diffusion 

(b) 



FIG. 6: (a) Schematic representation of a buckled grain 
boundary, and the induced periodic stress field alternating 
compression and traction regions in the adjacent crystals, (b) 
The grain boundary releases and captures crystal vacancies in 
the tractioned and compressed regions, respectively, produc- 
ing this way closed paths of the point defects. The diffusion of 
vacancies, equivalent to an atomic counterflow, have opposite 
senses in the two grains. 



As explained, the grain boundary elastic instability 
causes a strongly fluctuating stress field, which in the 
macroscopic scale averages to zero. Be 6b the energy 
necessary for the grain boundary to release a vacancy 
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into one of the two crystals adjacent to it, when the sys- 
tem is unstrained. When the grains are strongly stressed, 
this energy is expected to change linearly with stress, as 
long as e B is much larger than its expected variations. 
Denoting ajj the characteristic amplitude of the induced 
stress field spatial fluctuations, the relative concentra- 
tion of defects in two neighbouring regions of maximal 
and minimal stress is proportional to 



exp 



k B T 



+ P 1 ) 



exp 



k B T 

(40) 

where p 1 = — {a xx + a yy + a zz )/3 is the contribution from 
the induced stress field to the pressure. Expression (40) 
depends only on the normal stresses a xx , a yy and a zz 
because it is expected that the energy of a vacancy will 
not be significantly affected by the shear stresses, which 
do not change the volume of the hole, but only its shape. 

The flow of defects between these two regions is propor- 
tional to the spatial concentration variation ( 40 ) , which 
can be written as 



release and capture crystal vacancies, to or from the two 
grains, causing a relative motion of them. Also, in both 
theoretical approaches the in-plane shear force induces 
a strongly varying normal stress field localized in the 
boundary region and, as a consequence, the less ordered 
intergrain matter evaporates and condenses vacancies in 
a periodic series of sources and sinks. The main differ- 
ence is in the physical nature of the sources and sinks 
for the point defects. In Refs. and l2"2l the sources and 
sinks are originated by the buckling and trenching of the 
initially plane grain boundary (with eventual steps) by 
effect of the overcritical shear stress applied to it. Raj 
and Ashby |34j assume an uneven grain boundary sub- 
jected to a shear force. Sources and sinks are caused by 
the variation of the normal stresses induced in the oppo- 
site sides of the permanent grain boundary irregularities. 



V. THE CONSTITUTIVE EQUATION FOR 
SMALL STRAIN 



A. The equation 



knT 



e B + fl*p 



When the plastic distortion just starts, e ss and holds 

(41) 

the initial condition p8[), p = — a/3. Then 



for p 1 < p. In Eq. (41) D = D Q exp[-e a /(k B T)} is the 



diffusion coefficient for crystal vacancies, with e a being 
the corresponding activation energy. No matter the de- 
tails of the process, one can expect that the speed of the 
boundary sliding be proportional to the rate at which 
the assisting vacancies flow along the crystal boundaries, 
and expression (411 leads to the main Eq. (39 1 with 
£o = £s + e a - 

By the time, direct experimental evidence was supplied 
by Vetrano et al. [36 , 37 for the deformation induced su- 
persaturation of vacancies in the vicinity of many grain 
boundaries of an Al-Mg-Mn alloy during superplastic 
deformation. Rapid quenching of the sample while be- 
ing superplastically deformed revealed the formation of 
very unstable nano-cavities in a great fraction of the 
grain boundaries, which rapidly coalesce and disappear 
on moderate heating. The authors show they start form- 
ing after cooling and subsequent mounting of the sample 
in the transmission electron microscope, and conclude 
that the cavities provides evidence of a supersaturation 
of crystal vacancies near the sliding grain boundaries. 
The defects condense into voids when sliding abruptly 
ceases, and the just released vacancies are not recaptured 
by the grain boundaries. Hence grain boundary sliding 
is associated to the observed vacancy excess, which pro- 
vides evidence of the capital role of these defects in the 
deformation process. 

Some aspects of the mechanism for grain boundary 
sliding introduced in Refs . I2T1 and l22l and reviewed above 
in a more compact approach, resembles the one put for- 
ward earlier by Raj and Ashby [53]. In both the shared 
grain boundary is taken as a material medium able to 



1 . /2r c 
- arcsm - — - 

2 V a 



(42) 



and, replacing this and Eq. ( 39 1 in Eq. ( 36 ) 



--2C 



k R T 




2tc 
a 



x exp 



eo — r2*cr/3 
k~Wt 



(e«0). 



(43) 

Eq. ( 43 1 gives the plastic strain rate e of an axially 



symmetric sample at temperature T under the action of 
the uniaxial tensile stress a. Hence it is the constitutive 
equation of the material for the established conditions. 

The discussion on grain sliding of sections [H] and |III| 
is valid for any polycrystalline solid. Also, there is no 
reason to think that the grain boundary sliding mech- 
anisms leading to Eq. ( 39 1 for the coefficient Q be ex- 



clusive to superplastic grain sliding. Hence, in principle, 
Eq. ( 43 ) is expected to hold for the deformation of any 



fine grained polycrystalline solid, fine enough to ensure 
that grain boundary sliding be the rate limiting deforma- 
tion process. 

However, the large amount of experimental data avail- 
able in the literature on the superplastic materials makes 
them specially suited for testing Eq. (43 1. The data for 



very small strain rates, in the range 10 _o to 10 s , are 
specially valuable because the effect of the critical stress 
t c on the stress-strain rate curves is particularly notorius 
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in it. Such small strain rates are not usual in tensile tests 
of normal, not superplastic, materials. 



B. Comparison with previous results and with 
experiment 

To appreciate how t c affects the strain vs. strain rate 
curves, consider Eq. (43 1 in the asymptotic limit of very 
large tensile stress a, such that t c /<j — > 0. The equation 
reduces to 



krT 



eo — f2*<r/3 



(ewO, 



0), 



(44) 

which is the result given in Ref. [2TJ in which the effect of 
the critical stress r c discussed in sections |TT] and |fff| was 
passed over. This produces no significant consequence 
when comparing the theoretical results with data on al- 
loys having small r c , as is the case of the aluminium 
Al-7475 samples whose mechanical testing was reported 
by Hamilton [35]. However, the data on the same alloy 
Al-7475 prepared by Pilling and Ridley [39], aluminium 
Al-8090 SPF [3J or titanium Ti-6A1-4V [3J EHUD], do 
show some disagreement with Eq. ( |44[ ) for strain rates 
below lO^s" 1 . 

As an illustrative example, Fig. [7] shows the excellent 
fit given by Eq. p3| to the data of Cope et al. on Ti 
6A1-4V [31 EH 140] together with the curve representing 
Eq. ( |44| with the same values for the parameters. The 
latter has negative concavity over the whole range of i 
and does not exhibit the characteristic inflexion point, 
which is traditionally taken as a distinctive feature of 
superplasticity. However, both curves become practically 
the same for e > 5 x 10 _4 s _1 . 
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FIG. 7: Curves given by Eqs. (43 1 and (44 1, showing the 
effect of the critical shear stress t c . Values for the parameters 
are shown in the insets. Circles represent data of Cope et al, 
Refs. El [T0l gQl on the superplastic titanium alloy Ti-6A1-4V, 
strained at constant true strain rates and T — 1093 K. 

Figure [7] also illustrates the usefulness of very low 
strain rate data. As will be seen in the next sections, the 



critical shear stress r c has important consequences in the 
plastic deformation at finite strains and high or low strain 
rates, as fracture, for example, but does not produce any 
easily recognizable feature in the test curves. On the con- 
trary, at low enough e, tensile stresses are close to 2r c , 
and the solid angles A, B and C in Fig. [5] have signifi- 
cant magnitudes, producing the inflexion point shown in 
Fig. [7] As the deformation mechanisms considered here 
are in principle not exclusive to superplastic materials, 
we guess that the strain vs. strain rate curves for nor- 
mal ductile materials should exhibit and inflexion point 
as well, at low enough strain rates. 

In a subsequent paper [22], Eq. (44 1 was modified in 
the way 



n*(<x-<x ) 

C r-= ex P 

KrT 



knf 



(45) 



and an excellent agreement with the experimental results 
was achieved at low, intermediate and high i. However, 
the introduction of the threshold tensile stress <tq was 
based on just intuitive physical arguments, and justified 
only by the good results it produced. 




FIG. 8: The pre-exponential factor in Eq. (43 1 behaves like 
g — o"o with ctq = 2r c . 



Fig. [8] explains the success of Eq. (45 1, 
that, despite its more involved appearance, 



It shows 
the pre- 
exponential factor in Eq. (43 1 is almost linear in a, and 
vanishes for a = 2r c , which plays the role of <7q in 
Eq. (45 1. Therefore, Eq. (451 constitutes a very good 
approximation of the exact Eq. ( 43 1 , and the same ac- 



curate fit of the experimental data on Al-7475, Al-8090 
SPF [39], and Ti-6A1-4V E3 HQ] attained with the 
former [25] is to be expected for the latter, but now with 
no empirically justified modification. 

Fig. [9] shows the data of Cope et al. [3 EB HO] on the 
superplastic titanium alloy Ti-6A1-4V at six tempera- 
tures (circles). The solid lines represent Eq. (43 1. The 
constants fi*, eo and Cq have the same values for the 
six curves, and are shown in one of the insets. Notice 
that the values for f2*, eo and Co used here in relation to 
Eq. ( 43 1 are the same than those used in Ref. [22 to in- 



sert in Eq. (45 1, with the same success. The critical shear 
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stress r c evidences some dependence with temperature, 
as shown in the other inset. 
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FIG. 9: Stress-strain rate curves for Ti-6A1-4V at six tem- 
peratures. Circles represent data of Cope et ai, Refs. I51 llUI 
I40I and solid lines depict Eq. ( 43 1 for the parameters shown 
in the inset. The parameter r c exhibits some temperature 
dependence. 



where i is considered as a given constant. On the other 
hand, writting Eq. (36 1 in the form 



2di 1 



cot(20 c ) + 29 c 



7T 

2 



(47) 



and, derivating with respect to p with e constant, one 
has that 



2de 1 dQ ~, nn ,d0 c 



(48) 



Derivating also Eq. (39), the alternate expression for 
dQ/dp 



dQ 
dp 



n* 

knf' 



(49) 



is inferred. Replacing it in Eq. (48) and writting the 
result in differential form gives 



The critical shear stress r c is expected to be tempera- 
ture dependent. It may be equal to either the smaller of 
the two critical stresses associated to the grain boundary 
elastic instabilities studied in Ref s . |2"T1 and l2"2l or the crit- 
ical resolved shear stress (CRSS) of the material when it 
is greater than the critical stresses for the elastic insta- 
bilities. The critical stresses for the elastic instabilities 
depend on the elastic modulus of the material, which 
varies with temperature, and the CRSS is temperature 
dependent as well when the temperature is high enough. 
For Ti-6A1-4V it seems that r c is the CRSS. 



VI. FINITE DEFORMATION AT CONSTANT 
STRAIN RATE 



The strain function 



dp 
Q 



T ck B T o 



dfl*e 



T cot z {29 c )d9 c 



(50) 



Combining Eqs. (46 1 and (50 1, and integrating, one fi- 
nally obtains 



: = - M / d9 cot 2 {29) 



Bfl* 



ft, 



1 - cos(2(9) 



7rsin(26») 



2cos(26>) 7T 
n + 2 



sin(2(9) 
-i -l 



-26 



(51) 



where the limits 9 and 6 C correspond to the critical an- 
gles for strain e — and the final value e of the strain, 
respectively. 

Hence the strain e is related with the auxiliary variable 
9 C by an expression of the form 



Deformation up to a finite strain e can be accomplished 
in an indefinite number of inequivalent ways because of 
the other two variables involved in the process, a and 
e. Establishing the time dependence of one of them, or 
its dependence on the intermediate values assumed by e, 
would be enough to determine the process. The simplest 
conditions are to keep e — constant, as is usually done 
when testing superplastic materials, or set a = constant, 
which is the common practice for normal ductile materi- 
als. In this section we will assume the former. 



Eq. ( 34 1 can be put in differential form as 



de 



2de (I- cos(26» c ) 



29 

Bt c V sin(26» c ) 

46» c 2cos(26» c ) tt\ _1 dp 

7rsin(2(9 c ) 7t + 2J Q{p)' 



(46) 



(52) 



where F{9) is the universal function 



F{9) 



d9 cot (20) sin(20) 



(1 - cos(26>) 



29^{29)-^-^ + ^l 

IT 7T 2 



-i (53) 



whose graph is depicted in Fig. 10 F{9) is monotonically 
decreasing in its whole range (fJ, 7r/4), and diverges at 
9 = and 6 = 7t/4. 

The equation for the critical angle 9o at e = follows 
from combining Eqs. (42 1 and (43 1 in order to eliminate 
a. It is obtained 
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FIG. 10: The strain universal function F(9) and its asymp- 
totas at 9 — and 9 — it /A. 



B. The strain and strain rate dependent stresses 



Combining Eqs. ( 36 ) and ( 39 1 to eliminate Q, and then 



isolating p, one obtains the general equation 



k B T 

si* 



In 



( 2C n*r c 

V k B Te 
In fcot(20 c ) 



29 r 



SI*' 



(56) 



which gives the pressure once the relation 9 C = 9 c (s, e, T) 
is known from the procedure explained in the previous 
subsection. 

Knowing 8 C = 9 c (e, e, T) and p — p(e, e, T), the ten- 
sile stress a can be obtained from Eq. ( |37| , which can be 
written as 



2C 



SI*T C 



cot (26» ) + 26» 
1 

xexp -— — 
knT 



2Sl*~, 



Q) 



3sin(20 o ) 



(54) 



This equation gives 9q as a function of e and T. There- 
fore, the set of Eqs. ( [52] ) , (53) and (54 1 allows to deter- 
mine 9 C — 9 c (e, e, T). 



The adimcnsional coefficient appearing in Eq. ( 52 1 es- 
tablishes the scale for e and is a very small number. The 
bulk elastic modulus B is close to 0.7 x 10 11 Pa for alu- 
minium and 1.0 x 10 11 Pa for titanium at normal room 
temperature. From the comparison of our curves with ex- 
periments for e ss 0, we know that SI* is 2.6 x 10 -21 cm 3 
for Al-8090 and 5.9 x 10~ 22 cm 3 for Ti-6A1~4V at rather 
high temperatures, however, SI* is not expected to vary 
with T. Hence, at T = 300 K, 



k B T 
BSl* 



2.3 x 10~ 5 - 7.0 x 10~ 5 . 



(55) 



At elevated temperatures, B may decrease substantially, 
but anyway the coefficient ( 55 1 is expected to remain 
small. 



Because of the small value of the coefficient ( 55 1 , any 
significative strain e demands that either 9q or 9 C , or 
both, be in one of the two asymptotic regions 8 > 
or 9 < 7r/4. The asymptotic region 9 < 7r/4 can be 
associated to creep, while 8 > is related to ductile and 
superplastic deformation. For example, to attain a strain 
e ~ 1 - 10 % it must be F(9 C ) - F(9 ) ~ 1000. Observing 
Fig. [10] one readily realizes that such big differences can 
only exist in the asymptotic regions of F{6). 



4r r 



3sin(2(9 c ) 



(57) 



or, more explicitly, 



4r, 



k B T 



3sin(26l c ) SI* 



In 



2CqSI* t c 



k B Te 



•In 



(cot(20 c ) 



29 r 



_£p_ 
SI* 



(58) 



The transversal compression a x 
Eqs. (pi). 



<j y follows from 



C. The constitutive equation for a ^> t c 

The procedure outlined in the preceding subsection for 
determining the exact dependence of the flow stress on 
strain, strain rate and temperature relies on solving a 
set of trascendental equations. However, a rather simple 
closed-form constitutive equation a = cr(s,s,T) can be 
written if one finds acceptable to restrict the scope to 
the strain rates for which Eq. (44 1 is valid. This means 



e > 2 x 10 s _1 for Ti-6A1~4V, which is a typical figure 
of more or less general validity. 

The coefficient k B T/(BSl*) appearing in Eq. (52 1 is 
in general a very small adimensional quantity. As F(6) 
is an universal function exhibiting large variations only 
close to its sigularities at 8 — and 9 — tt/4, significant 
strains e demand 9 C > or 8 C < 7r/4. In the latter 
eventuality the stress is near the minimum necessary to 
produce plastic flow, which progresses at a very low rate 
and corresponds to creep. The former case is the ductile 
deformation we are interested in. 

When 9 r > the universal function F(9), given by 



Eq. ( 53 1 , can be substituted by the asymptotic equation 



F{9) 



TT — 8/lT l 



(59) 
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and then the expression ( 52 I for the strain becomes 



k B T 1 



8/tt Bn* e c 



(60) 



where £o is an integration constant determined by the 
condition ( 38 1 . Isolating 8 C from Eq. ( 60 ) and inserting 



in Eqs. (56 1 and (57 1, both rewritten for 9 C <C 1, it gives 



k B T 

n* 



In 



(n-8/Tr)C (n*) 2 BT c 
{k B T) 2 e 



(e + £o) 
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(61) 



and 



2(7r-8/7r) BO*r c , 
g= 3 fc.T C ( g + £ o)-P- (62) 



The initial condition (38), that is, p = — cto/3 when e = 0, 
with cr being the flow stress for e — 0, allows us to put 
Eqs. dwj and ([62) as 
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In 



C ft* f(n -8/ir)Bn*T c 
k B Te 



k R T 



Cq 



" ft* ' 

(64) 

which is the general constitutive equation a — cr(s, i, T) 
we were searching for. If one replaces in Eq. (64 1 e = 



and a = ao, the approximate constitutive equation (44) 



for e — is recovered. The domain in which the approxi- 
mation holds is apparent in Fig. 7] where the curves given 
by Eq. ( 44 1 and the exact Eq. ( 43 1 are compared. 



D. Exact expression for the slope da/de at e = 

Derivating and combining properly the equations writ- 
ten in subsections lVI Al and lVI Bl one can deduce a closed- 
form equation for the slope da/de at e — 0. In general, 
for any e, 



da 

de 



=B 



An* 



i 



3fc B Tcos(26» c ) cot(20 c ) + 20 c - tt/2 
1 - cos(26» c 



sin(26i c ) 
2cos(26» c ) 7T 

7T + 2 



29 A 1 



7rsin(26> c ) 



(65) 



For e = one can replace sin(26> c ) = 2t c /<7, and Eq. (65) 
becomes an expression of only the initial value of a. 



The slope da/de at e = is a very important quantity 
because it is observed that always the maximal super- 
plastic strain to failure is attained close to a temperature 
and strain rate such that 



9a 

de 



0. 



(66) 



which may be considered as an analytical condition for 
superplasticity. Although this conclusion is mostly em- 
pirical, it is quite reasonable because, as a polycrystal 
cannot flow steadily by the monotonic increase of the in- 
ternal pressure p, Eq. (66 1 expresses the situation which 
best emulate a steady regime. The next example illus- 
trates the point. 
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FIG. 11: The uniaxial flow stress a for a modified Al 5083 
alloy deformed at 570 °C and constant strain rate e. Circles 
represent the data of Kaibyshev et. al. (Ref. I41|> and solid 
line depicts Eq. ( 43 1 with the material constants shown in 
the insets. 



Fig. [TT] shows the variation of flow stress with strain 
rate for a modified 5083 aluminium alloy, deformed at 
T = 843 K and constant true strain rate. The experi- 
mental points are due to Kaibyshev et. al. [41] and the 
solid line depicts Eq. (43 1 with the constants shown in 
the insets. These authors observed a sharp maximum in 
elongation to failure of 1150% for T = 843 K and a strain 



rate e = 2.8 x lO^s -1 . 

Inserting in Eq. (65 1 the constants obtained from the 



fit of the stress-strain rate data, shown in Fig. [TT] the 
slope da/de of the stress-strain curves at e = can 
be calculated for any initial stress a. The solid-dashed 
straight lines drawn in Fig. [12] represent the calculated 
tangents to the stress-strain curves at e = start- 
ing from a = 1.33, 2.4, and 7.0 MPa. According to 
Eq. (43 1 these three stresses correspond to the strain 
rates e = 2.8 x 10~ 3 , 1.4 x 10" 2 and 7.0 x 10 -2 s _1 . 
The circles in Fig. [T2] represent the experimental stress- 
strain data taken by Kaibyshev et. al. at constant strain 



rates e = 2.8 x 10^, 1.4 x 10~ 2 and 7.0 x 10~ 2 s~ 



The 



consistency between theory and experiment is apparent 
in Fig. 12 The effective bulk modulus, which fixes the 
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FIG. 12: Circles represent stress-strain data at constant 
strain rates i = 2.8 x 10~ 3 , 1.4 x 10~ 2 and 7.0 x 10~ 2 s _1 
for the same system of Fig. |11| The straight lines are the tan- 
gents to the stress-strain curves at e — and the same strain 
rates, calculated with Eq. (651. The effective bulk modulus is 
B = 3000 MPa. 



scale for e, was taken as B = 3000 MPa. Best superplas- 
tic response of the material is attained at a strain rate 
e = 2.8 x 10~ 3 , for which da/de « at e = 0. 



E. Grain sliding and the sign of the concavity of 
the stress— strain curve 



In a recently published experimental study of the 
plastic tensile deformation of Ti-6A1-4V, Vanderhasten 
et. al. find out that positive concavity seems to be the 
fingerprint of grain boundary sliding 42 , 43 . These au- 
thors conducted tests starting from room temperature 
to 1323 K, recording tensile properties, microstructural 
evolution and crystallographic texture changes at a con- 
stant true strain rate e = 5x 10 -4 s -1 . The tests seem 
not intended to attain optimal superplastic performance, 
but clearly show a dramatic increase of the strain to 
fracture at T = 873 K (600° C), which reaches the op- 
timum (e = 1.6) at T = 1123 K (850° C). Starting from 
T = 998 K, grain boundary sliding begins to be dominant 
and, coincidently, the corresponding stress-strain curve 
exhibits positive concavity over the most of the strain do- 
main. The lower temperature range is characterized by 
negative concavity and grain stretching. At temperatures 
between T = 1123K and T = 1173 K the stress-strain 
curves have positive concavity over the whole range, from 



e = to fracture at e ~ 1.6. For higher temperatures the 
strain to fracture decreases and the stress-strain curves 
develop a final sector with small negative curvature be- 
fore failure. 
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FIG. 13: Circles represent data for Ti-6A1-4V of Vander- 
hasten et. al. at the the temperatures for which grain sliding 
dominates ove r grain stretching. Solid lines were obtained 
from Eq. ( |64[ ) and the same set of parameters, which are 
shown in the insets. Just Co shows some variation with T, 
probably due to its grain size dependence. The effective bulk 
modulus is B = 3000 MPa. 
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FIG. 14: Stress-strain rate curves calculated from Eq. (431 
and the same material constants of Fig. |13| No experimental 
data are available to compare with. 



The open circles in Fig. 13 represent experimental data 
taken by Vanderhasten et. al. on the stress-strain re- 
sponse of the titanium alloy Ti-6A1-4V at the three tem- 
peratures for which the deformation is observed to take 
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place by grain sliding 



The strain rate was kept at 



e = 5.0 x 10 s in all runs. The samples exhibit high 
ductility but not properly superplasticity. The solid lines 
represent the asymptotic expression (64 1 for the three 
temperatures, with the parameters shown in the insets 
and B = 3000 MPa. Eq. ( 64 I can describe with high pre- 



cision the experimental situation, which can be shown by 
comparing the values for a at e = given by Eq. ( 64 I with 
the exact value ctq, obtained from Eq. (43) to introduce 
it into the same Eq. (64) as a one of its constants. The 



agreement between theory and experiment is remarkable 
up to e = 1. 

Fig. 14 shows the stress-strain rate curves given by 
Eq. (43 1 and the same material constants used to fit 
the data of Fig. |13| There is no experimental points 
to compare with since Vanderhasten et. al. employed a 
unique strain rate. The curves exhibit the typical sig- 
moidal shape. 



Ti-6AI-4V 



n"= 1.73 X 10" 21 [cm 3 ] 
1.85 [eV] 



T c =0.35 [MPa] 




° T[K] CqIs 1 ] q 

1123 8800 -0.237 
1173 10400 -0.237 
1223 8000 0.000 



0.5 



1.0 
STRAIN S 



1.5 



may be as follows: The mean grain sizes measured by 
Vanderhasten et. al. at 1123 and 1173 K show a rather 
irregular reduction on strain. The consequent division 
of old grains involves the creation of new grain surfaces 
which do not have necessarily the composition and struc- 
ture of the old ones. Grain boundaries are efficient traps 
for impurities, and can condense atomic species having a 
low bulk concentration. However, the process is governed 
by bulk diffusion and is not rapid. Therefore, concerning 
their composition and consequent sliding ability, fresh 
grains may be very different than well-aged ones. From 
this viewpoint, grain size variations may simply reduce 
the number of grains able to slide. This explains also 
why the smooth stress-strain curves do not reflect the 
irregularities observed in the grain size-strain curves. 

The theoretical approach to finite strain superplastic- 
ity put forward here permits to explain quantitatively 
the main features observed in the superplastic deforma- 
tion. However, the theory lacks an important element: 
the effect of heat generation at the sliding grain bound- 
aries. In our scheme, the whole work made by the ex- 
ternal forces is dissipated into heat at the sliding grain 
boundaries, which occupy a very small volume when com- 
pared with the grain bulk. Hence one may expect that 
grain boundaries will be at a higher temperature than 
bulk matter, being the nominal temperature the bulk 
one. Grain boundaries whose normals are within solid 
angles, determined by 9 Cl around the principal stress di- 
rections do not slide because the in-plane shear stress is 
smaller than r c , and then remain at the bulk tempera- 
ture. Bending downward of the stress-strain curves at 
high strains and strain rates seems to be a grain bound- 
ary temperature effect. 



Two phase alloys and the nature of the 
superplastic deformation 



FIG. 15: Sa me as Fig. 13 but with Co replaced by Co, as 
given by Eq. (67 1 with the values of q shown in the inset. Cq 



models the effect of grain size variations. 



In general, Co depends on the grain size d as d~ 3 
and the rest of the parameters do not vary with d |22j . 
Fig. [15] shows an attempt to explain the departure be- 
tween theory and experiment, starting at e = 1, by the 
dynamic modification of the grain size. The curves rep- 
resent Eq. ( 64 1 with Co substituted by 



C'o 



(1 + ge 3 )C , 



(67) 



where q is a constant. The values of q giving the best fit 
to experiment are shown in one of the insets of Fig. |15| 
Amazingly, q turns out to be such that C = at the 
fracture strain. This means that failure occurs because 
the grains become unable to slide, which can hardly be 
explained by just their size variations. The explanation 



The scheme described thus far should hold for sin- 
gle phase alloys. Materials combining two or more al- 
lotropic phases in thermodynamic equilibrium may dis- 
play a more complex phenomenology. In particular, su- 
perplastic materials are generally two-phase or pseudo 
single phase alloys. Probably the most characteristic ex- 
ample of the former is Ti-6A1-4V, which combines an 
a hep solid phase with a (3 phase with bec structure, in 
comparable volume fractions up to the (3 transus tem- 
perature. Cope el al. [3l [10] observed best superplastic 
response of Ti-6A1-4V for a temperature T — 1153K, 
for which the (3 phase proportion was 0.42%. In pseudo 
single phase materials the second phase occurs as small 
inclusions between the grains constituting the main one. 

The new feature introduced by the coexistence of two 
allotropic phases in thermodynamical equilibrium is that 
a volume reduction is not necessarily related with just 
the elastic properties of the material. If the two phases 
(which we will call generically as a and (3) have different 
specific volumes, the total volume may change also by the 
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transformation of part of one of the phases into the other 
one. In general, the atomic fraction of, for instance, the 
(3 phase 



N 



(68) 



where Np is the number of atoms of the main constituent 
of the alloy in the phase and N = N a + Np is the 
total, is a funtion of the temperature T and pressure 
p. If p changes with time, the atomic fraction of the 
phases changes, as well as the total volume V. Hence, 
the dilation rate is no more given by Eq. ( 33 1 but by 



V 
V 



1 



V{3 



B N a v a + Npv/3 \ dp 



ON, 



P, 



(69) 



where v a and vp are the volume per atom of the two 
phases. 



Eq. ( 69 1 assumes that the atomic fractions have their 
equilibrium values in any instant, and hence a quasi- 
static process, otherwise the transformation kinetics 
would come into the fore. The total area of the interfaces 
between the phases is maximal when their volumes have 
the same value, and then comparable volume fractions 
contribute to increase the phase transformation rates, 
and the accuracy of Eq. ( 69 1 . Also, the new isother- 



mal coefficient appearing in Eq. ( 69 ) is in general not a 



constant, and it is expected to be much more sensitive to 
temperature than 1/B. With this in mind we define 



The abnormally large strains to failure displayed by 
superplastic materials have been attributed to the better 
stability they have against necking in tensile tests. The 
plastic flow can be represented by the relation 



a = Ke r ' 



(72) 



where m is a function of i and e, generally known as the 
strain rate sensitivity, which is always smaller than unity. 
When ttl — 1 any irregularity in the cross section of the 
sample is not accentuated by the deformation in a tensile 
test, necking has no effect and the material behaves as 
a Newtonian viscous fluid. Necking instability is high 
for m below 0.5 and decreases markedly when m > 0.5. 
Indeed, the condition of high values of m is a necessary 
requirement, but not enough to reach elongations as large 
as those termed superplastic. Our previous discussions 
show that plastic elongation has a limit because demands 
a volume reduction, and failure is not a matter of just a 
mechanical instability of geometric origin. 

Therefore, in our view, superplasticity has its real ori- 
gin in three concurrent features of the material and phys- 
ical conditions of the deformation process: (a) A high 
strain sensitivity m. (b) The ability of the material to re- 
duce the internal pressure p by undergoing a phase trans- 
formation to a phase of higher density, (c) The material 
must be in the verge of a ductile to brittle transition. 
This latter requirement is explained below, taking the 
superplastic behaviour of titanium Ti-6A1-4V to illus- 
trate the argument. 



1 



and call 



Vf) ~ Va 
N a V a + Npvp 



1 1 1 

B* B B a p 



dp 



(70) 



(71) 



B a n is a positive quantity since the two factors in the 



right hand side of Eq. (70 1 have opposite signs, as de- 



manded by mechanical stability. 

The equations of the previous sections are formally re- 
covered this way, with an effective bulk modulus B* in- 
stead of the purely elastic one B. 

The elastic modulus B takes in general large values 
(typically B ss 1 x 10 5 MPa for metals), which is not nec- 
essarily the case of B a p. For two phases coexisting in 
thermodynamical equilibrium one may expect that the 
volume fraction of each one would vary significantly for 
pressure variations much smaller than the yield stress. 
Then, values of less than 1 MPa for B a p are quite rea- 
sonable. In such a situation 1/B <C 1/B a p and 1/B* 
is several orders of magnitude larger than 1/B. By this 
reason solids exhibiting two allotropic phases in equilib- 
rium may display a special plastic behaviour, of which 
superplasticity seems to be an example. 



G. The high temperature second yield point 

The general properties of the stress vs. strain curves 
can be visualized from just examining Eq. ( 64 1 for the 



strain dependent flow stress cr. It has a first term linear 
in e, with slope proportional to 1/T, and a negative log- 
arithmic term whose coefficient increses with T. Hence, 
at low temperatures cr is linear in e and, as T increases, 
the slope decreases and the function develops positive 
concavity. Beyond a critical temperature for which the 
slope at the origin vanishes the curve starts at e = 
whith negative slope and, after surpassing a minimum, 
or second yield point, begins to increase monotonically. 

Hence, Eq. (64 1 predicts strain hardening at any strain 
for temperatures below a critical temperature, which 
roughly coincides with the temperature which optimizes 
superplastic properties. At higher temperatures the ma- 
terial should weaken for strains between e = and the 
second yield point. This configures a mechanically unsta- 
ble situation which cannot occur in practice in a smooth 
and homogeneous way. The equations written so far de- 
scribe the mechanical response at a point, which does 
not necessarily can be straightforwardly extended to the 
whole sample. Homogeneity is a quite reasonable as- 
sumption to study averaged properties in the macro- 
scopic scale, but may fail in critical situations. We can 
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safely assume that f2* and eo are uniform and constant 
throughout the material, but Co depends on the grain 
size as dr^ [22 , and hence may fluctuate or have a mul- 
timodal distribution, following grain size. 

There is a class of mild steels [44 and aluminium al- 
loys [45j |46] exhibiting stress-strain curves with two yield 
points. The plastic deformation of these materials is in- 
homogeneous and takes place by a series of succesive 
sudden band formations (Liiders bands [E]), transver- 
sal to the tensile axis and making a rather precise an- 
gle with it (of about 56°), which concentrate the strain. 
The flow stress keeps constant (though displaying a ser- 
rated aspect) until the bands, increasing in number and 
extension, fill the whole sample and the complete ma- 
terial becomes strained beyond the second yield point 
strain. Just then the stress starts to increase. In the 
intermediate steps, samples have sections with no strain 
and sections strained beyond the second yield strain. The 
relevant question is then how a system with the stress- 
strain curves given by Eq. (64) will evolve with strain, 



and what kind of effective stress-strain curves could be 
expected from such evolution. 

If our material is not exactly homogeneous the critical 
temperature for the onset of the second yield point is not 
the same in its whole extension. Hence, at temperatures 
close to these critical temperatures and small strain one 
may expect the solid will be partitioned in a network 



of sectors hardening on strain, and another one which 
weakens. As the stress is uniform, the weakening sectors 
will deform much faster, and will be responsible for the 
main part of the effective strain before reaching their 
second yield point. Meanwhile the hardening sectors will 
furnish mechanical stability against sudden yield steps, 
like Liiders band generation. 

Anyway, the condition for superplastic behaviour 
seems related with the onset of the second yield point, 
producing a transition to a potentially unstable system. 
Such a transition begins at the temperature for which 
da/de at e — vanishes and changes from positive to 
negative sign. On the other hand, generally speaking, 
the most basic condition for a steady flow, as superplas- 
tic deformation tends to be, is that the dynamical vari- 
ables remain constant and independent of e, which here 
represents time. As a varies with e, all what can be done 
is to emulate a steady flow by providing the conditions 
that keep da/de small and steady in a range as large as 
possible, which incorporates e = 0. 



Acknowledgments 

The authors thanks the support by the Programa Bi- 
ccntcnario de Ciencia y Tecnologia (Chile) ACT-26. 



J. Pilling and N. Ridley, Superplasticity in crystalline [17 
solids, (The Institute of Metals, Camelot, Southampton, 
UK, 1989). [18 
T. G. Nieh, J. Wadsworth and O. D. Sherby, Superplas- 
ticity in metals and ceramics (Cambrige, UK 1997). [19 
M. T. Cope and N. Ridley, Mat. Sci. Technol. 2, 140 
(1986). [20 
T. H. Thomsen, D. L. Holt and W. A. Backofen, Met. [21 
Eng. Quart. 2, 1 (1970). [22 
Aerospace Materials, edited by B. Cantor, H. Assender [23 
and P. Grant. (Institute of Physics Publishing, Bristol, [24 
UK, 2001). 

M. A. Khaleel, H. M. Zbib and E. A. Nyberg, Int. J. [25 
Plast. 17, 277 (2001). [26 
M. B. Taylor, H. M. Zbib and M. A. Khaleel, Int. J. Plast. [27 
18, 415 (2002). 

D. W. Livesey and N. Ridley, Metall. Trans. 9A, 519 
(1978). [28 
D. W. Livesey and N. Ridley, Metall. Trans. 13A, 1619 [29 
(1982). [30 
M. T. Cope, D. R. Evetts and N. Ridley, J. Mat. Sci. 21, 
4003 (1986). [31 
A. Ball and M. M. Hutchinson, Met. Sci. J. 3, 1,(1969). 
A. K. Mukherjee, Mater. Sci. Eng. 8, 83 (1971). [32 
T. G. Langdon, Phil. Mag., 22A, 689 (1970). 
H. W. Hayden, S. Floreen and P. D. Goodall, Metall. [33 
Trans. 3A, 833 (1972). 

M. F. Ashby and R. A. Verrall, Acta Metall. 21, 149 [34 
(1973). [35 
R. C. Gifkins, Metall. Trans. 7A, 1225 (1976). 



A. Arieli and A. K. Mukherjee, Mater. Sci. Eng. 45, 61 
(1980). 

H. Fukuyo, H. C. Tsai, T. Oyama and O. D. Sherby, ISIJ 
International 31, 76 (1991). 

C. G. Wang, M. W. Fu, C. X. Cao and H. B. Dong, Mat. 
Sci. Eng. A, in press. 

M. Lagos and H. Duque, Int. J. Plast. 17, 369 (2001). 

M. Lagos, Phys. Rev. Lett. 85, 2332 (2000). 

M. Lagos, Phys. Rev. B 71, 224117 (2005). 

M. Lagos, Phys. Rev. B 73, 224107 (2006). 

P. Bellon and R. S. Averback, Phys. Rev. Lett. 74, 1819 

(1995). 

D. J. Srolovitz, Acta Metall. 37, 621 (1989). 

J. Muller and M. Grant, Phys. Rev. Lett. 82, 1736 (1999). 
M. A. Grinfeld, Dokl. Acad. Nauk. SSSR 265, 836 (1982); 
Sov. Phys. Dokl. 31, 831 (1986); Europhys. Lett. 22, 723 
(1993). 

D. R. M. Williams, Phys. Rev. Lett. 75, 453 (1995). 

F. Y. Genin, J. Appl. Phys. 77, 5130 (1995). 

D. E. Jones, J. P. Pelz, Y. Hong, E. Bauer and I. S. T. 

Tsong, Phys. Rev. Lett. 77, 330 (1996). 

D. E. Jesson, K. M. Chen, S. J. Pennycook, T. Thundat 

and R. J. Warmack, Phys. Rev. Lett. 77, 1330 (1996). 

P. C. Searson, R. Li and K. Sieradzki, Phys. Rev. Lett. 

74, 1395 (1995). 

R. Q. Hwang, J. C. Hamilton, J. L. Stevens and S. M. 

Foiles, Phys. Rev. Lett. 75, 4242 (1995). 

R. Raj and M. F. Ashby, Metall. Trans. 2, 1113 (1971). 

Y. J. Wei and L. Anand, J. Mech. Phys. Sol. 52, 2587 

(2004). 



18 



[36] J. S. Vetrano, E. P. Simonen and S. M. Bruemmer, Acta 
Mater. 47, 4125 (1999). 

[37] J. S. Vetrano, C. H. Henager and E. P. Simonen, in Super- 
plasticity - Current Status and Future Potential, edited 
by P. B. Berbon, M. Z. Berbon, T. Sakuma and T. Lang- 
don. (Materials Research Society, Warrendale, PA, 2000). 

[38] C. H. Hamilton, C. C. Bampton and N. E. Paton in Su- 
perplastic Forming in Structural Alloys, edited by N. E. 
Paton and C. H. Hamilton. (The Metallurgical Society of 
AIME, Warrendale, PA, 1982). 

[39] J. Pilling and N. Ridley, in Aluminium Technology '86, 
edited by T. Sheppard. (Institute of Metals, London, UK, 
1986) 

[40] M. T. Cope, M. Sci. Thesis, Victoria University of 
Manchester, U. K. (1982). 



[41] R. Kaibyshev, F. Musin, D. R. Lesuer and T. G. Nieh, 

Mater. Sci. and Eng. A 342, 169 (2003). 
[42] M. Vanderhasten, L. Rabet and B. Verlinden, J. Mat. 

Eng. Performance 16, 208 (2007). 
[43] M. Vanderhasten, L. Rabet and B. Verlinden, Materials 

and Design 29, 1090 (2008). 
[44] R. Song, D. Ponge and D. Raabe, Acta Mater. 53, 4881 

(2005). 

[45] J. S. Hayes, R. Keyte and P. B. Prangnell, Mater. Sci. 

and Technol. 16, 1259 (2000). 
[46] Z. C. Wang and P. B. Prangnell, Mater. Sci. and Eng. A 

328, 87 (2002). 
[47] R. Hutanu, L. Clapham and R. B. Rogge, Acta Mater. 

53, 3517 (2005). 



